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TECHNICAL PAPER 


MECHANICS OF LIQUID HELIUM IN A PARTIALLY FILLED ROTATING DEWAR 
IN LOW GRAVITY - WITH APPLICATION TO GRAVITY PROBE-B 

INTRODUCTION 


The Gravity Probe-B (GP-B) spacecraft (Fig. 1) is designed to test the Equivalence Principle 
through long-term (1 year) monitoring of forces on a set of gyros in free-fall around the Earth. 
Extraneous forces on these gyros must be kept at very low levels, corresponding to accelerations of 

lO* 1 ^ g or less. This will require a drag-free (to 10‘ 10 g) control system, which uses a proof mass 
similar to the experiment gyros as its sensing element. In addition, the measurement system for moni- 
toring these low force levels will require cooling to below superconducting transition temperatures. 

The approaches to both cooling and control involve the use of liquid helium. The boil-off from the 
liquid helium dewar will be used as a propellant to maintain the drag-free operation of the spacecraft. 
The requirement for an operational lifetime approaching one year means that a large quantity of liquid 
helium must be used, and that it will be gradually depleted over the experiment lifetime. This varying 
amount of liquid helium gives rise to the possibility of several problems which can degrade the GP-B 
experiment. Most of these potential problems are due to asymmetry in the static liquid helium dis- 
tribution or to waves in the free surface which will be present in the partially-filled dewar. Before 
discussing these problems in detail, however, it must be determined whether classical fluid mechanics 
applies to the analyses. This determination is necessary because of the variety of unusual effects 
observed in liquid helium at low temperatures. 


FLUID MECHANICS OF LIQUID HELIUM 

Helium exists in nature as a mixture of two isotopes, He^ and He^. The relative abundance of 
the He ^ is only about one part in 1 0^, however [ 1 ] . Helium can then for most purposes be considered 
to consist of He^ only. This isotope, when in a liquid state, has two phases called Hel and Hell. 

Figure 2 [2] is the low temperature phase diagram for He^ (helium). There is a phase transition line, 
called the X-line, which separates the liquid Hel and liquid Hell phases. The terminus of this line on the 
gas-liquid line is called the X-point. This denotes the temperature at which the Hel-Hell transition 
occurs at saturated vapor pressure conditions. This temperature is 2.172°K. 

Liquid Hel is found to behave normally. That is, experiments on liquid Hel exhibit properties 
and phenomena which are common to viscous fluids. The theory which describes this behavior is 
classical fluid mechanics. Liquid Hell (i.e., liquid helium below 2.172°K) is found experimentally to 
exhibit behavior which cannot be explained classically. Some of the properties observed are: 

1 ) The fountain effect. This effect is the flow of liquid helium through capillaries (or through 
pores) from colder to hotter regions. 




Figure 1 . Gravity Probe-B. 



Figure 2. Low temperature phase diagram for He^. 
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2) The mechano-caloric effect. This is the inverse of the fountain effect. Here, when liquid 
helium flows through porous material from one reservoir to another, the reservoir which is losing 
material experiences a rise in temperature. 

3) Ability to flow through narrow regions such as slits and capillaries without sustaining a drop 
in pressure. This effect disappears above some critical opening size and some critical velocity. The 
opening size and velocity critical values are dependent on each other. 

4) Related to this last effect, in experiments involving the rotation of liquid helium in 
geometries below critical size and at subcritical velocities, the liquid helium is seen to only partially 
contribute to the moment of inertia of the container. At the X-point, all the liquid helium participates 
in the rotation, but at lower temperatures, the measured moment of inertia of the container and liquid 
helium decreases as a function of temperature. Interpretation of these results requires a two-fluid model 
of liquid helium involving a normal viscous component and a superfluid component. 

For the GP-B spacecraft liquid helium management problem, the geometries are so large that the 
critical values will be exceeded. The liquid helium is treated, then, as a normal viscous fluid (with a 
very low viscosity, on the order of 15 micropoise). This has been experimentally verified for large 
containers by Mason, et al. [3]. The laws and methods of classical fluid mechanics will therefore apply 
here. 


GP-B EXPERIMENT CONFIGURATION 


The planned configuration for GP-B is shown in Figure 1. The region of the spacecraft dewar 
containing the liquid helium is nearly an annular volume between concentric cylinders. In discussing 
liquid helium configuration and dynamics that is the model which will be used (Fig. 3). The helium 
liquid and gas will be assumed to occupy the annular region of length 294 cm, outer radius 54 cm, and 
inner radius 18 cm. 



HELIUM 

(LIQUID AND GASEOUS) 


INNER CYLINDER 




r 2 = 54 cm 


= 18 cm 



Figure 3. GP-B model system. 
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LIQUID HELIUM CONFIGURATION 


The requirement that accelerations experienced by the proof mass be kept at the 10 "*^ g level 
or lower means that the liquid helium used for cooling and propulsion can be a source of problems for 
the GP-B system. The proof mass responds to external forces as the sensing portion of the drag com- 
pensation system. When the proof mass experiences external forces, the gaseous helium propulsion 

system responds with thrusts which will null net forces to the 10"*® g level or lower. If the forces on 
the proof mass are due to spacecraft drag, then the desired effect is achieved. If the forces are from 
other sources such as gravitation, then the effects of accelerating the spacecraft can still null net forces 
on the proof mass, but not on the experimental gyros since the effects of the gravitational forces will 
vary over the spacecraft. Also the effect of accelerating the reference frame (the spacecraft) to com- 
pensate for a force internal to the frame will result in orbital changes for the spacecraft. 

In order for the GP-B drag compensation system to function properly, internal forces must be 
kept below the 10"*® g design requirement. 

Gravitational forces can perturb the proof mass if there are distributions of liquid helium which 
are not both axially symmetric and longitudinally symmetric about a plane containing the center of 
mass and perpendicular to the rotation axis. Liquid helium distributions lacking these symmetries will 
result in gravitational effects on the proof mass in two different ways. First, the unsymmetric liquid 
helium distribution will have a mutual attraction with the proof mass. If we consider worst case mass 
distributions, we can establish upper bounds on these gravitational attractions. For the dewar half full, 
and with all the liquid helium at one end, the gravitational acceleration produced at the proof mass is 

of the order of 10'^ g or the dewar half full, with the liquid helium on one side, the resulting accelera- 
tion level at the proof mass is about 10"^ g (see Appendix A). 

The second problem introduced by unsymmetric liquid helium mass distributions is the gravity 
gradient effect. Distributions of liquid helium which lack the proper symmetry (as described above) 
will cause the center of mass of the spacecraft to shift away from the center of the proof mass. The 
path of the proof mass in orbit will then be slightly different than that followed by the spacecraft 
center of mass. The force on the proof mass resulting from this nonfree-fall trajectory can be computed 
by expanding the expression for the Earth’s gravitational attraction about the spacecraft center of mass. 
After doing this, we find that the magnitude of accelerations due to this effect at the proof mass can be 

of the order of 10"^ g or greater for worst case mass distributions (see Appendix B). 

We have established that unsymmetric distributions of liquid helium can result in acceleration 
levels at the proof mass which can be more than two orders of magnitude above experiment design 
requirements. Whether or not these distributions can arise depends on the balance of forces acting on 
the liquid helium. Forces which must be considered are: 

1) Spacecraft drag 

2) Spacecraft propulsion 

3) Gravity gradient forces 

4) Gravitational forces between the liquid helium and the spacecraft 
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5) Centrifugal forces due to spacecraft rotation 

6) Surface tension forces in the liquid helium. 

For a spacecraft operating as designed, spacecraft drag will be balanced by the helium propulsion 

system. There should be no net accelerations (above the 10’*® g level) on the liquid helium from these 
sources. In Appendix C, the expected magnitudes of gravitational, gravity gradient, centrifugal, and 
surface tension forces on the liquid helium are compared. Gravity and gravity gradient forces are shown 
to be at least one order of magnitude smaller than the smallest expected centrifugal forces. For 
purposes of discussing the liquid helium configuration then, spacecraft drag and propulsion and gravity 
and gravity gradient forces can be neglected. To a close approximation the liquid helium configuration 
will be the result of a balance between centrifugal forces and surface tension for the boundary condi- 
tions imposed by dewar shape. 

When the GP-B spacecraft is deployed it will not be spinning. In the early stages of the experi- 
ment, a spin rate of up to about 0.1 rad/sec will be imposed for instrument calibration. After calibra- 
tion, the rotation rate will be reduced to its operational value of approximately 0.01 rad /sec. Changes 
in rotation rate of the spacecraft are not felt immediately in the liquid helium. Angular momentum 
must be diffused and convected into the liquid helium from the side and ends of the dewar. Until the 
helium is in solid body rotation with the spacecraft, pressure distributions due to centrifugal forces will 
deviate from hydrostatic and will be unpredictable. Appendix D deals with this problem. It is shown 
there that spin-up times are very short (a day) compared to experiment lifetime. 

Knowing that equilibrium configurations of the liquid helium are governed by the centrifugal 
force-surface tension balance, and that spin-up can be accomplished fairly rapidly, permits analyses 
which can make some predictions about expected configurations. These analyses use two approaches, 
energy analysis and pressure balance. 

Energy analyses are used to predict whether, for a given set of conditions, a simply connected 
(“spherical”) bubble of gaseous helium or a multiply connected (“toroidal”) bubble is expected. 
Appendices E and F deal with this approach. Results show that for the operational rotation rate, and 
for high fill factors (small bubble volumes), conditions are marginal for forming the desired toroidal 
bubble. That is, an unsymmetric liquid helium distribution could result from these conditions. 

Pressure balance calculations have been performed in order to obtain equilibrium liquid helium 
distributions as a function of rotation rate and fill factor. Appendix G uses pressure balance to obtain 
gaseous helium bubble shapes in a longitudinal section containing the rotation axis. Here curvature is 
considered in the r-z plane. Appendix H assumes that an axisymmetric toroidal bubble exists, and 
computes equilibrium configurations as a function of rotation rate. In the calculations in Appendix 
H, the helium gas-liquid interface is assumed to have a finite radius of curvature only in the r-9 plane 
(i.e., in a plane perpendicular to the rotation axis). Results from these pressure balance computations 
confirm that axisymmetric liquid helium configurations are not expected for high fill factors at the 
operational rotation rate of 0.01 rad/sec. 


CONCLUSIONS 


Liquid helium mass distributions in the GP-B dewar are possible which are capable of producing 
accelerations on the proof mass on the order of 10'^ g. Since the proof mass is displaced from the 
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experiment masses, these accelerations cannot be compensated by the drag-free propulsion system. Also, the 
response of the drag-free propulsion system to these kinds of internal forces could degrade the orbital 
dynamics of the spacecraft. It was shown that the configuration of the liquid helium in the dewar is 
primarily a function of centrifugal forces and surface tension. Energy and force balance analyses 
indicate that at the planned operational rotation rate, the desired axisymmetric liquid helium configura- 
tion is improbable for high fill factors and marginal for lower fill factors. Approaches to achieving the 
symmetric distribution will be discussed in a separate report. 
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APPENDIX A 


WORST CASE MASS DISTRIBUTIONS AND THEIR EFFECTS 


Worst Case Longitudinal Displacement of GP-B Center of Mass 

The worst case GP-B center of mass displacement due to an asymmetrical configuration of the 
liquid helium will exist if the dewar is half full and all the liquid helium is at one end (Fig. A-l). 


GASEOUS HELIUM 

LIQUID HELIUM 


GASEOUS HELIUM 

LIQUID HELIUM 


Figure A-l. Worst case — axial. 

Assuming a total dry mass Mj of the dewar to be 1921.8 kg, the mass M 2 of liquid helium 

(for completely filled dewar) to be 320.5 kg, and the dewar length to be 294 cm, the center of mass 
displacement along the axis for this worst case is given by: 


Mj(0) M 2 (294/4) 
c.m. Mj + m 2 Mj + M 2 


(160.25) (73.5) 
1921.8 + 160.25 


= 5.66 cm . 


Worst Case Radial Displacement of GP-B Center of Mass 

The worst case of center of mass displacement in a radial direction for the GP-B will occur if 
the dewar is half full and all the liquid helium is at one side (Fig. A-2). 
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V 


’I" ■ H, _ ' 


Figure A-2. Worst case — radial. 
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Assuming a uniform mass density, p, of liquid helium, the location of the helium center of 
mass will be given by: 




7T 

= - COS 8 | 

0 




= 4_ (r2 3 - r l 3 > 
3n ( r 2 2 ' r l 2 ) 

for rj = 18 cm and ^ = 53 cm 



= 24.43 cm . 


This value of the liquid helium center of mass displacement gives a radial displacement of the 
total GP-B center of mass of: 
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_ Mj(0) M 2 (Y c m ) 
AR c.m. M j + M 2 + Mj + M 2 


(160.25) (24.43) 
1921.8 + 160.25 


c.m. = 1.88 cm . 


Gravitational Attraction Between the Proof Mass and the Liquid Helium 

The displacement of the liquid helium center of mass from the spacecraft center of mass (hence 
from the proof mass center of mass) results in gravitational attraction between the liquid helium and the 
proof mass. This can be expressed as an acceleration at the proof mass location due to the liquid 
helium: 


a = 


GM 

T 


where M is the liquid helium mass and r is the distance between the proof mass and liquid helium 
centers of mass. 

The maximum center of mass displacements for a half full dewar have been calculated. For the 
longitudinal case this results in the following acceleration level: 


GM (6.67 x 1 O' 8 ) (1.60 x 10 5 ) , 7 

a.j = = cm/sec^ 

r 2 (7.35 x 10 1 ) 2 


= 1.98 x 10' 6 cm/sec 2 = 2.02 x 10' 9 g . 


For the case where the dewar is half full with all the liquid helium at one side: 


(6.67 x 10' 8 ) (1.60 x 10 5 ) , 7 

an = cm/sec“ 

(2.443 x 10 1 ) 2 


= 1.79 x 10' 5 cm/sec 2 = 1.83 x 10' 8 g . 
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APPENDIX B 


GRAVITY GRADIENT EFFECTS 


The acceleration due to the gravitational attraction of the Earth at a point on an orbiting space- 
craft may be expressed as an expansion about the spacecraft center of mass : 



where G is the gravitational constant, M is the mass of the Earth, and R is a vector from the center of 
the Earth to a point on the spacecraft. Now let 


R = Rj + R2 - (R e + h) Rj + R2 * 


where Rj is the vector from the center of the Earth to the spacecraft center of mass, R2 is the vector 

from the spacecraft center of mass to the point in question, R e is the radius of the Earth, and h is the 
spacecraft altitude (Fig. B-l). 



Figure B-l. Gravity gradient. 
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Ill I 


Also, 


^ _ R _ R 1 + R 2 
R 

|R| |R| 

We can then write: 


= - GM 


R 1 + r 2 

l R l 3 


Rl + R2 

[R • R] 3 / 2 


or 


3g = - GM (Rj + R 2 ) [(Rj + R 2 ) • ( R i + r 2 )]' 3/2 
= - GM (Rj + R 2 ) [Rj 2 + 2 Rj • R 2 + R 2 2 ]* 3 / 2 


= - GM (Rj + R 2 ) 


Ri 


/ I+2 V!i + !Hij 

\ Rj 2 Rj 2 / 


-3/2 


= - GM (Rj + R 2 ) 


, R 1 3 


1 + 2 


R 1 ' r 2 R 2" 


-3/2 


R 


1 Ri 


This can be expanded in a binomial series. If we wish to have accuracy to the third power of 
1 /Rj , we need to retain only the first two terms in the expansion. We have then: 


a g = - GM 



✓x 

R 


- 3 


1 ‘ R 2 
Ri 
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Again keeping terms through 1/Rj , 


a g a - GM 


/N 

Ri 


r 2 2 


/\ — 

R, • R 


- 3 


1 • *2 £ 


R- 


R l 3 


Ri + 


Ri 
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We observe that the first term is just the gravitational acceleration at the spacecraft center of 
mass. This term is balanced by the centrifugal acceleration of the spacecraft. For a spacecraft in an 
inertial flight mode, the second and third terms represent unbalanced accelerations. The result will be 
that a mass at any point away from the spacecraft center of mass experiences gravity gradient forces. 
The component of displacement perpendicular to the radius generates forces directed back toward the 
center of mass. The component of displacement along the radius rj generates forces away from the 
center of mass. 

For GP-B, then, if there is any displacement of the spacecraft center of mass from the proof 
mass center of mass, gravity gradient forces will act on the proof mass in the manner described above. 
Any distribution of liquid helium lacking axial symmetry (including symmetry about a plane perpen- 
dicular to the axis and through the proof mass) will cause these forces. It is useful to make worst 
case estimates of the magnitude of these effects. 

If R -2 is taken to be a vector to a point fixed in the GP-B spacecraft, and is assumed to lie in 
the orbital plane, then the angle between Rj and can be represented by cot + 0 where co is the 

orbital angular rate, t is the time, and is a fixed place angle. If we take 4 > = 0 for the angle at which 
Rj and R 2 are aligned, and let t = 0 at that point, then we have: 


/\ _ _ 

Rj * R_2 = IR2 1 cos cot 


The gravity gradient part of the acceleration can then be written : 


a = - GM 


r 2 /\ r 2 

- 3 cos cot Ri + 

Rj 3 1 R } 3 


= - GM 



r 2 a 

cos cot Ri 

V ' 


+ r 2 


when R 2 is parallel to Rj : 


?(0) ■ - GM — - [- 2 Rj] 


where R 2 is antiparallel to Rj : 
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R 2 /v 

a(7r) = - GM — - [2 Rj] • 

R 1 

When R 2 is perpendicular to Rj, i.e., cos cot = 0, we have: 

r 2 

a(7r/2) = a(37r/2) = - GM — - [R 2 ] . 

R 1 3 

Examination of these expressions for gravity gradient acceleration shows that it is directed 
toward the center of mass of the spacecraft for "rj i r^, and away from the center of mass for rj 

parallel to r 2 , and has maximum magnitudes of 2 GMR 2 /R 1 3 . Where: 

G = 6.67 x 10" 8 cm 3 gm'^ sec' 2 

M = mass of Earth = 5.98 x 10 27 gm 
Rj = Re + h 

R e = radius of Earth = 6.37 x 10 8 cm 

h = altitude = 520 km = 5.2 x 10 7 cm 
a = gravity gradient acceleration 

= 2(6.67 x 10' 8 ) (5.98 x 10 27 ) (R 2 )/(6.89 x 10 8 ) 3 
= 2.44 x 10' 6 R e sec' 2 = 2.49 R 2 x 10' 9 g , 

where R 2 must be expressed in cm, and represents the amount of center of mass displacement due to 
asymmetries in the helium configuration. 
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APPENDIX C 


FORCES ON LIQUID HELIUM 


Gravity Gradient Accelerations 

From Appendix B, the expression for the gravity gradient acceleration at a point represented by 
the vector R2 from the spacecraft center of mass is given by: 



An estimate of the magnitude of this acceleration was shown to be : 


a = 2.5 R 2 x 10" 9 g , 


where R 2 is expressed in centimeters. For R 2 = 150 cm, 


a = 4 x 10' 7 g . 


Gravitational Accelerations 

The gravitational acceleration everywhere inside the GP-B dewar can be determined in principle 
as : 


"3(r) = - V 0(f) ’ 


where 0(r) is the gravitational potential: 


-/■ 


,/ . . P(r ) 3 , 

0(r) = / d J r 

|r - r‘| 


Here p(x5 is the mass distribution of the spacecraft. 

To determine whether it is necessary to evaluate this integral, we may obtain an estimate of 
a(r). An upper bound on the magnitude of a(r) will be given by assuming all the mass of the spacecraft 
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is concentrated at the center of mass and computing the acceleration at a distance L/2 away, where L is 
the dewar length. For L = 150 cm, and M (the spacecraft mass) to be 1900 kg, we have : 


a 


G M 
~ 


_ (6.67 x 10'^ cm^ gm"* sec" 2 ) (1.9 x 10^ gm) 
(1.5 x 10 2 ) 2 

= 5.63 x 10'^ cm/sec 2 
= 5.74 x 10" 9 g . 


Surface Tension Forces 

The pressure difference across any fluid-fluid interface (Fig. C-l) is given by: 



where a is the interfacial tension measured as a force per unit length, and Rj and R 2 are the principal 
radii of curvature at a point on the interface. 



Figure C-l. Pressure across 2-fluid interface. 

For liquid helium various values for a have been reported. A value on the upper end of the 
range is 0.53 dynes/cm [1]. Other values range downward to about 0.12 [3]. Using 0 = 0.53 dynes/ 
cm and a characteristic radius of 18 cm (the inner cylinder radius for GP-B), we obtain AP = 0.029 

dynes/cm . A pressure can be calculated which would result from a hydrostatic head due to an acceler- 
ation field of 4 x 10" 2 g (the value derived as an upper limit of gravity gradient effects). This would be 
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P = p g h = (0.145 gm/cm 3 ) (4 x 10' 4 cm/sec^) (54 cm) 


= 0.0031 dynes/cm 7 , 


where the radius of the GP-B dewar provides the 54 cm length scale. 

This value is a factor of 10 smaller than the surface tension forces. We can expect then that 
surface tension will be much more important than gravity gradient effects in determining liquid helium 
configuration. 


Centrifugal Forces 

Since the GP-B spacecraft will rotate about its longitudinal axis, the liquid helium in the dewar 
will experience centrifugal forces. The acceleration field due to this rotation is given by : 

2 

a c = orr , 

where co is the rotation rate. At the operational rotation rate of about 0.01 rad/sec, we have : 
a c = 10‘ 4 r . 

Expressed in effegravity this becomes : 
a c = 10' 7 r g , 
where r is in centimeters. 

For r = 18, a c = 1.8 x 10'^ g and for r = 54, a c = 5.4 x 10"^ g. These levels are about an 

order of magnitude higher than the largest expected gravitational effects and comparable to surface 
tension effects. 
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APPENDIX D 


SPIN-UP AND SPIN-DOWN OF THE GP-B 


A fluid in a container undergoing a uniform rotation about a single axis will ultimately come 
into equilibrium with the container and the container-fluid system will rotate as a solid body. When 
this condition is achieved, a hydrostatic pressure distribution exists in the fluid. The pressure distribu- 
tion may readily be computed from the centrifugal force and is just a function of distance from the 
rotation axis, fluid density, and rotation rate. When the fluid is not spun-up, that is when it is adjusting 
to a changing angular velocity of the container, the pressure distribution is very complex. In order to 
have an experimental condition which is known, it is preferable for the liquid helium in the GP-B 
dewar to be spun-up while the experiment is being conducted. It is necessary then to compute 
estimates of the time required for the liquid helium to adjust to changing angular velocities. This char- 
acteristic time will be denoted as the spin-up time. 

The model appropriate for the GP-B dewar is concentric cylinders with the liquid helium 
contained in the annular region between them (Fig. D-l). The problem then is to get angular momen- 
tum (vorticity) into the liquid helium from the rotating cylinders. There are two transport processes 
available to carry this vorticity generated at the container walls (see W. B. Watkins, R. G. Hussey, Spin- 
up from rest in a cylinder, Physics of Fluids, Vol. 20, No. 10, Pt. 1, Oct. 1977). First, vorticity gen- 
erated at the walls of the dewar diffuses into the liquid helium. The characteristic time for this process 

is T( j ~ a 2 /V, where a is the radius of the outer cylinder and v is the kinematic viscosity of the liquid 
helium. For a = 53.34 cm and for v = ju/p = dynamic viscosity/density = (1.3 x 10'^ poise)/0.1452 
gm/cm 2 = 8.95 x 10‘ 2 cm 2 /sec, we have a characteristic time: 

= (53.34 cm) 2 /8.95 x 10' 2 cm 2 /sec 

= (28.45 x 10 2 )/8.95 x 10’ 5 sec 

= 3.18 x 10 2 sec = 368 days « 1 year . 



Figure D-l. Spin-up model. 
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Spin-up through diffusion of viscosity from the walls is then a very slow process, requiring a time of the 
order of the experiment lifetime to occur. The second process which contributes to spin-up is the con- 
vective transport of vorticity due to centrifugal pumping (Eckman pumping) at the cylinder end walls. 
Figure D-2 indicates the nature of this effect. Liquid in contact with the end wall must rotate at the 
angular rate, w, of the end wall. Shear stresses in the viscous liquid entrains a layer having thickness 

the order of 5= O/w)^ ( see fj P. Greenspan, The Theory of Rotating Fluids, Cambridge Univ. 

Press, 1968, p. 32). Here 5 = (8.95 x 10"^ cm^/sec)/10"^ rad/sec ^ = 9 5 x 10"^ C m, that is, 6 is 
about 1 millimeter. The fluid in this layer has angular velocity and therefore is subject to centrifugal 
forces which push it radially outward. The sidewalls turn the flow. For a closed cylinder where there 
is centrifugal pumping at each end, toroidal flow cells develop, extending from each end to the midplane 
of the cylinder. The spin-up time scale due to this mechanism is given by: 


r c = L(vc c)" 1 ^ 


where L is the cylinder length. For L - 294 cm this gives: 


t c = (294 cm) (8.95 x 10'^ cm^/sec) (10'^ rad/sec)"^/^ 
= 3.12 x 10~* sec = 3.61 days . 



Since the spin-up time is directly proportional to the length scale L, a decrease in this time may 
be obtained by making L smaller. This could be done by using baffles which would be plates parallel 
to the cylinder ends. Then the characteristic length would be the baffle spacing. If this spacing is 
about L/10, then the characteristic spin-up time reduces to about 1/3 day. 

Since r c « Tj, the spin up will be completely dominated by centrifugal pumping. r c is short 
enough that hydrostatic conditions should be established quickly (compared to total experiment time). 

Spin-down occurs in a fluid by the same mechanisms as spin-up. Therefore, reductions in rota- 
tion rate for GP-B will be accomplished through the end wall pumping effects, and will have the same 
characteristic time as for spin-up. 
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APPENDIX E 


ENERGY ANALYSIS 


Given a vapor bubble of specified volume, we will eventually be concerned with its position both 
along the axis and in the radial direction. However, for the time being, ignore the bubble’s position 
along the axis and consider only its radial position. For the case where the cylinder is not rotating, the 
bubble might look like Figure E-l. As long as the bubble was small enough to fit between the 
cylinders, its shape would be spherical. (Note that liquid helium at this temperature is a perfectly 
wetting liquid thus preventing the bubble from attaching and forming a solid-vapor interface.) With zero 
rotation, the radial position of this vapor bubble would be arbitrary, but with even a slight rotation the 
bubble would begin to float toward the center of the container as the heavier liquid was forced toward 
the outside (Fig. E-2). 



Figure E-l. Bubble in nonrotating case. 



Figure E-2. Bubble with small co. 

When the bubble reaches the inner cylinder its final shape depends on a combination of factors. 
Those parameters which will affect the bubble’s shape and bubble size for a given container include the 
density difference between the vapor and liquid phases, the vapor pressure, the surface tension of the 
vapor liquid interface, and the rotation rate. When the rotation rate is small relative to the surface 
tension, the bubble would look something like Figure E-3. However, as the rotation rate increases (or 
the surface tension decreases) the bubble will deform, wrapping itself around the cylinder as in Figure 
E-4. This deformation would continue with increased rotation until the two ends of the bubble touch. 
At that time the bubble would “pop” and reform into a toroidal shape (Fig. E-5). 
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Figure E-3. Bubble with c o' increasing. 




Figure E-4. Bubble, large to. 



Figure E-5. Axisymmetric bubble. 

It is easy to imagine that if we then slowed the rotation rate, the vapor would remain in some 
torus-like form, perhaps even until the cylinder stopped rotating completely. The significance is that 
for one rotation rate there may be two very different possible states, one axisymmetric and one not. 
Which shape is thermodynamically preferred depends on the balance between the surface energy and 
the potential energy due to rotation. The surface energy would be minimized by a spherical bubble 
off to one side while the potential energy due to rotation would be minimized by a bubble wrapped 
around the inner cylinder. Unfortunately, it is extremely difficult to calculate the actual shapes of these 
two states and to subsequently compare their energy. Instead we have chosen to simplify the problem 
by comparing the energy difference between the undistorted shapes of a sphere and a torus. The 
question then is what is the energy required to go from a torus to a sphere (Fig. E-6). 
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AE 


Figure E-6. Energy model. 


The change in energy in going from a torus to a sphere of the same volume is the sum of the 
energy released by the decrease in surface area plus the energy required to transport the liquid helium 
through the acceleration field to its new position. (The energy required to relocate the vapor is 
neglected.) The first part of this sum is easily calculated and is simply a(A g -Aj) where a is the surface 
tension, A s is the surface area of the sphere, and Aj is the surface area of the torus. The final result 
is that the difference in energy between a torus and a sphere is: 


n 

AEtotaL = ff ( A T ” A s ) P ft 
(Torus -+ Sphere) 


(W« 1+ b)[R lt 2 b J-(j|^ + Rl |.4) 


where 


a = surface tension 
Aj = area of the torus 
A = surface area of the sphere 
co = rotation rate (radians/sec) 
p = density (g/cm^) 
b = small radius of the torus 
Rj = radius of the inner cy Under 
a = radius of the sphere . 

If AE is positive for a given bubble volume and co then the bubble would rather be a sphere. If E 
is negative then a torus is the preferred configuration. A value of AE = 0 implies that the bubble is 
at a crossover point where it would just as soon be a torus as a sphere. 
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In order to get a feel for what shape might be expected, some examples have been worked for 
helium in the dewar with the dimensions shown in Figure E-7. 



Figure E-7. Energy model with GP-B parameters. 

Using the values of a = 0.37 erg , p = 0.145 g/crn^, Rj = 17.8 cm, and w = 0.01 rad/sec, a 

bubble with a volume of 24 liters (1/100 of the available volume of the actual dewar) would rather be 
spherical (AE = 975 dynes/cm). However, an increase to a rotation above 0.38 rad/sec would mean that 
the bubble would prefer to be a torus (AE < 0). Similarly, a vapor bubble with a volume of 510 liters 
at a rotational speed of 0.01 rad/sec would rather be toroidal, but if slowed down to below W = 0.0089 
would prefer to be a sphere except that such a large sphere could not fit undistorted in the space 
between the inner and outer cylinder. 

The conclusion of this energy analysis is that near a rotation speed of 0.01 rad/sec (the pro- 
jected rotation rate for the actual experiment) it is questionable what form the bubble will take. 
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APPENDIX F 


SURFACE ENERGY AND ROTATIONAL POTENTIAL ENERGY 


A slightly different energy argument from that of Appendix G is presented here. Two possible 
configurations of the liquid helium will be analyzed with respect to surface energy and rotational poten- 
tial energy. These configurations are an axisymmetric LHe distribution (Fig. F-l) and a simply con- 
nected bubble (Fig. F-2) of helium gas. The simply connected bubble may or may not extend the 
length of the dewar. 



Figure F-l. Axisymmetric bubble. 



Figure F-2. Simply connected bubble. 

The configuration with the lowest possible rotational potential energy is just that shown in 
Figure F-l. In this configuration, all the potential energy has been converted to rotational kinetic 
energy. Ignoring the contribution of the helium gas (which is negligible due to its very low density of 

the order of 10' J times the LHe density), the rotational kinetic energy of this configuration is 


,27T 


E R = f f f \ P T 2 dV 
e=0 z=0 r 2 


1 

4 


7 r 




co 2 (r 3 4 - r 2 4 ) 
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where £ is the length of the dewar. To obtain the rotational potential energy of any configuration, we 
must find the difference between the rotational kinetic energy of that state and E^. 

The total energy of a given configuration will be the sum of its rotational potential energy and 
its surface energy. Surface energy is just : 

Eg = crA > 

where A is the total interfacial area. The total energy is then: 


E T E s + ^potential 
= a A. + Ej^ - Ej^' , 

where E^' is the rotational kinetic energy of the given state. 

For an axisymmetric bubble, the surface energy is: 

E s c = a 2 tt r 2 £ . 

The total energy for this configuration, Ej c is then: 

E T C = E s C + % ~ • 

Here, however, % ~ e r' so that for the axisymmetric (cylindrical) configuration: 
E-p c = 2ir X2 £ o . 


For the simply connected bubble (which will be labeled “spherical”) shown in Figure F-2, the 
rotational kinetic energy is: 


A6 £ , 2 7r £ r, 

E * S = / / / 2" r2 “ 2dV+ / /„ / - 2 ^ 2 “ 2 ^ ■ 


0 hr. 


A6 0 r, 
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Evaluating this integral gives: 


E r s = ip co 2 (r 3 4 - q 4 ) [2 tt£ - A0h] , 

where h is the length of the gas bubble and Ad is its angular extent. 

The rotational potential energy of this configuration is then: 

^potential “ E r " E r s 

= X - it £ p a; 2 (r 3 4 - r 2 4 ) - E r s 

= ^ 7r £ p o? (rj 4 - r 2 4 ) + ^ Ad h p co 2 (r 3 4 - rj 4 ) 

The surface energy of this configuration is : 

E s s = ctA = a ^2(r 2 - rj) h + ^ a A d (r 3 2 - rj 2 ) > 

where 

a = 1 for bubble stuck to end 
a = 2 for bubble away from end . 

The total energy for the simply connected bubble is: 

e t s = V + Er - E r s 

= a j^2(r 2 - rj) h + ^ a Ad (r 3 2 - rj 2 )J- ^ n £ p co 2 (r 2 4 - q 4 ) + ^ A0 h p co 2 (r 3 4 - rj 4 ). 

The criterion for having an axisymmetric cylindrical bubble is that the energy for that state be lower 
than that for alternate states: 

e t c < e t s 
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or 


E t s - E t c > 0 . 

This condition is: 

a |^2(r 2 - rj) h + i a Ad (r 3 2 - q 2 ) - 2i r r 2 

+ i Ad h p to 2 (r 3 ^ - r^) > 0 . 

8 

This is subject to the condition that the volume of t] 
y c = v s 
or 

IT £ (r 2 2 - rj 2 ) = i (A6 h) (r 3 2 - rj 2 ) 

so that: 

2t r £ (r 2 2 - q 2 ) 

(Ad h) = . 

(r 3 2 -ri 2 ) 

Assume that h is some fraction of £: 

h = b £ where 0 < b < 1 

then 

2tt (r 2 2 - rj 2 ) 

Ad = — . 

b (r 3 2 - rj 2 ) 


Rewriting the energy condition we have: 
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a j^2(r 2 - rj) h + ir (r 2 2 - q 2 ) - 2?r r 2 £ J - ^ ir ft p co 2 (r 2 4 - q 4 ) 


+ i tt £ p co 2 (r 2 2 - rj 2 ) (r 3 2 + rj 2 ) > 0 


or 


F = a ^2(r 2 - rj) h 


+ (J) " (r 2 2 -ri 2 )- 2 7T r 2 £ 


+ i 7T £ p CO 2 (r 2 2 - I-] 2 ) (r 3 2 - r 2 2 ) > 0 


Then, F(r 2 , a, b, co) > 0 defines the region in parameter space where an axisymmetric bubble 
is energetically preferred. Figures F-3 and F-4 show these regions for the following parameters: 

rj = 18 cm 

r 2 = variable (set for each calculation) 
r 3 = 54 cm 
£ = 294 cm 

b = a constant between 0 and 1 
a = 1 or 2 

co = independent variable 
p = 0.145 gm/cnr* 
a = 0.53 dynes/cm . 
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ENERGY BALANCE CRITERION FOR VARIOUS VALUES OF ULLAGE 
VOLUME, V (BUBBLE AT END OF DEWAR) 


Figure F-3. Energy condition, bubble at end of dewar. 



ENERGY BALANCE CRITERION FOR VARIOUS VALUES OF ULLAGE 
VOLUME, V (BUBBLE AWAY FROM END WALLS) 


Figure F-4. Energy condition, bubble away from end of dewar. 




APPENDIX G 


PLOTTING THE BUBBLE SHAPE 


The curvature of a bubble and therefore its shape is dictated by its surface tension and the 
pressure difference between the inside of the bubble and the outside. The relevant equation is: 


AP = - o (Kj ± K 2 ) , 


where AP is the pressure difference between the inside and outside of the bubble, a is the surface 
tension and Kj and K 2 are the orthogonal curvatures (Fig. G-l). Consider a nonrotating container of 
liquid helium containing a bubble of volume V (Fig. G-2). The pressure of the vapor inside the bubble 
is a function of the temperature. As heat flows into the container and vaporizes the liquid helium the 
overall temperature is held constant by means of evaporation through a porous plug. Since the tem- 
perature is fixed the vapor pressure inside the bubble likewise remains at a fixed value. The pressure 
in the liquid outside the bubble is then given by : 


^Liquid “ a ^1 ± ^2^ + ^vapor • 

For a sphere Kj ± K 2 = 2/R where R is the radius of the sphere. Therefore for a sphere : 

^liquid “ ^vapor " 2cr/R . 

Since the vapor pressure is a constant, the pressure of the liquid must increase over the course of the 
experiment as helium is vaporized and vented overboard (Fig. G-3). 



er - — « (k, s k 2 > 


Figure G-l. Pressure drop across bubble wall. 
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Figure G-2. Bubble in closed container. 


30 




Pi 


_ P 2 

0 

p,<p 2 

O 


Figure G-3. Pressure change with bubble growth. 

For a rotating system the problem becomes more complex but the equation Pjjq - P ya p = 
-a (Kj ± K 2 ) still holds. For the configuration given in Figure G-2 the Navier Stokes equations in 
cylindrical coordinates can be applied to get an expression for the pressure in the liquid: 


P 


3V r 3V r V ff 3V_ 3V r 

+ V r + + V 9 — 

3t r 3r r 30 2 3z 




d\ t 1 dv, t 1 
g f 2 r 9r r 2 g#2 


v r 2 iv e 

* dz 2 " r 2 “r 2 * e 


For the case where the fluid rotates as a solid body around the axis this equation reduces to : 



where p is the density, the velocity in the theta direction, P the pressure, and r the radial coordi- 
nate. Substituting cor = (where co = radians/sec) and integrating gives: 

2 2 
p co^ r z 

Pressure liquid = + Constant . 

2 


Substituting this expression into P^q - P yap 


Curvature (Kj ± K^) = - 


2 2 
p co^ x L 

2 ct 


= - a (Kj ± K 2 ) results in: 
Constant + Py a p 0r 

* 

o 
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For a “two-dimensional bubble” K 2 = 0. The resulting bubble can be plotted using a computer 
plotting program. An example of such a bubble plot is given in Figure G-4. The figure gives an 
example for two different values of the constant in the above equation. The curvature for these two- 
dimensional bubbles with K 2 set equal to zero is sharper than would be the curvature of the cross- 
section of an actual bubble. 

For the three-dimensional bubble which is axisymmetric (toroidal in shape) the value of 1/r 
can be substituted for K 2 again simplifying the problem where it can be easily graphed by computer. 

(This is an approximation for K 2 which will be exact at the inner cylinder and at the outer extremity of 

the bubble. The approximation should only affect the detailed shape of the bubble ends, not its size.) 
Figure G-5 is the side view of a torus drawn by the computer program. In the figure, S is the surface 
tension; o> is the rotation rate; the Int Press (8000) is the pressure of the helium vapor; and Ext Press 
is the pressure of the liquid helium. The bubble shown in Figure G-5 would grow as the helium was 
depleted and at some point would look like Figure G-6. The axial position of the bubble in Figure 
G-6 is for the most part arbitrary. The center of gravity of such a bubble is unpredictable and 
therefore undesirable. One possible remedy to this problem would be to step up the rotation rate. 

A higher rotation rate is more likely to produce an axisymmetric torus and will also stretch the bubble 
out such that at a certain volume the bubble would be positioned by the end walls. Figure G-7 shows 
the growth of a bubble at a higher rotation. The bubble shown in Figure G-7 would touch both end 
walls when its volume reached approximately one-fourth of the total volume. 



Figure G-4. Bubble shape, rd section. 
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Figure G-5. Bubble shape, rz section. 



RUN #1 S - .37 DYNES/cm W - .01 RAD/SEC. D«. 145 9/on INT PRESS - 8000 (DYNES/an 2 ). 


Figure G-6. Bubble shape, co = 0.01 rad/sec. 
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Figure G-7. Bubble shape, co = 0.05 rad/sec. 








APPENDIX H 


JL 


CYLINDRICAL BUBBLE PRESSURE BALANCE 


The pressure balance at a cylindrical interface, concentric about the rotational axis of the liquid, 
for the system shown in Figure H-l is given by: 

o ^ + 0 ^ = ^ (P2 " Pi) w2 ( r2 ” r l 2 ) 

or 

= (P2 ~ Pi) <^ 2 r2 ~ (P2 " Pp <^ 2 r l 2 


= (P 2 - Pi) u 2 r3 - (P 2 ~ Pi) ^ 2 r l 2 R 


2 a 
~R 

2 a 


R- 5 - r, z R - 


2a 


= 0 


(p 2 ~ Pi) ^ 



Figure H-l. Cylindrical bubble model. 

This can be solved numerically to obtain R as a function of angular rate co. R will be the radius 
of a cylindrical bubble which is in equilibrium at a given rotation rate. Results are shown in Figure H-2 
for values of surface tension of 0.35 and 0.53 dynes/cm. 

Another way to present the pressure calculations is to compute pressure drop across an inter- 
face as a function of radius of curvature of a toroidal bubble and to compare this to the pressure due 
to centrifugal forces at the same radius. For surface tension calculations are: 


AP = a/r2 • 


For centrifugal forces, calculations are: 
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AP = “ (p 2 - Pl> <^ 2 O^ 2 - iq 2 ) 


where P 2 = density of liquid helium and pj = density of gaseous helium. These calculations are plotted 
in Figures H-3 through H-5. 

It is necessary to determine how full the liquid helium dewar can be at a given rotation rate and 
still have an axisymmetric bubble. This can be done by first finding the bubble radius consistent with 
equilibrium (Fig. H-2) and then determining what fill factor this represents (fill factor is 1-c, where c is 
the ullage fraction). Figure H-6 presents plots of ullage fraction versus bubble radius. Curve 1 of Figure 
H-6 represents the GP-B system as designed. Curve 2 illustrates how an axisymmetric bubble could be 
maintained for larger fill factors if a larger diameter dewar is used. Figure H-7 shows the ullage fraction 
which can be maintained at a given rotation rate assuming cylindrical bubbles. 



ROTATION RATE (RAD/SEC) 


Figure H-2. Equilibrium cylindrical bubble radius versus rotation rate. 
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Pressure (dynes/c^) 



Figure H-3. Pressure versus radius of curvature for cylindrical bubble. 
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Figure H-4. Pressure versus radius of curvature for cylindrical bubble. 
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